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Abstract 

Advanced, multidimensional models are typically applied when researching processes occurring in the nearshore. Relatively 
simple, empirical equations are commonly used in coastal engineering practice in order to estimate extreme wave run-up 
on beaches and coastal structures. However, they were mostly calibrated to the characteristics of oceanic coasts, which have 
different wave regime than a semi-enclosed basin like the Baltic Sea. In this paper we apply the formulas to the Polish Baltic 
Sea coast. The equations were adjusted to match local conditions in two test sites in IVli^dzyzdroje and Dziwnowek, where 
beaches are under continuous video surveillance. Data from WAM wave model and coastal gauge stations were used, as 
well as precise measurements of the beaches' cross-sections. More than 600 run-up events spanning from June to 
December 2013 were analysed, including surges causing dune erosion. Extreme wave run-up /?2% was calculated and 
presented as a percentage value indicating what part of the beach was inundated. The method had a root-mean-square 
error of 6.1 and 6.5 percentage points depending on the test site. We consider it is a fast and computationally undemanding 
alternative to morphodynamic models. It will constitute a part of the SatBahyk Operating System-Shores, delivering forecasts 
of wave run-up on the beaches for the entire Polish coastline. 
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Introduction 

Coastal zone represents a dynamic environment where land, sea 
and atmosphere meet. It is also a system where up to two-thirds of 
the world population lives and many other are visiting those areas 
frequently [1,2]. In order to sustainably manage recreational, 
industrial, agricultural and any other activities taking place near 
the shores it is necessary to understand the dynamic nature of both 
marine and terrestrial interactions at the coasts. 

Small-scale elements of a depositional coast can experience 
rapid change in their morphology, mainly due to short-term sea 
level variations caused by storm surges [3-6] . In this perspective it 
is crucial for coastal managers to determine the impact of the sea 
on the coast especially in the aspect of potential hazard and risk to 
human activities as well as to the shore itself. This includes also the 
possibility of making predictions of those impacts. 

Recently, an evaluation of risks, existing coastal management 
plans and civil engineering-produced protection schemes in nine 
European countries revealed that operational approaches to 
coastal risks and hazards had been virtually non-existent in most 
of those countries [7] . The biggest effort used to be channelled 
only to forecasting physical parameters (wind, waves, sea levels 
etc.), without considering their impact on the coast. One of the first 
steps to adopt an alternate approach to coastal hazards was the 



MICORE project in which nine prototype early warning systems 
based on open-source morphological XBeach model [8] were 
developed [9,10]. MICORE made it possible to test, with good 
results, the XBeach model at a variety of test sites along the 
European coasts, including the Polish southern Baltic coast [11]. 

Advanced, multidimensional models such as XBeach are 
typically applied when researching processes occurring in the 
nearshore. They mathematically explain typically complex and 
non-linear relationships between different variables and are 
indispensable in enhancing knowledge about coastal systems 
[12,13]. Models provide support in finding the best means to deal 
with as yet unforeseen events and their effects, both morphological 
and socio-economic. In Poland, efforts are being made to expand 
research on other parts of the coast. Currently three test sites are 
included in the SatBahyk Operating System-Shores [14], the 
coastal-research component of a research project dedicated to the 
monitoring of the Baltic Sea [15]. Unfortunately, these models 
have a major drawback: they are computationally demanding, 
limiting the possibility to use them for detailed analysis of large 
parts of the shoreline. 

But alongside those robust models some simpler methods were 
developed. At times when solving complex, non-linear equations 
was problematic to say the least, coastal engineers created 
empirical formulas. They were intended to improve the designs 
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of coast protection structures. Wave run-up could be calculated 
and therefore enabled the engineers to design proper heights of 
embankments or seawalls [16,17]. In the 1980s scientists started 
field experiments trying to apply the equations to natural beaches, 
an attempt that proofed successful. Several new erjuations were 
also developed [18-24]. Yet unlike more sophisticated methods, 
they are not being used in operational forecasting, but are rather 
confined to engineering practice. 

Nevertheless, the reliability of aforementioned methods, both 
the empirical equations and models, needs always to be verified 
when applied to a new area. One possibility is to conduct in situ 
measurements of the desired phenomena. Much more convenient, 
though, is to use an increasingly popular alternative in the form of 
video monitoring systems. This method was applied for coastal 
research in the early 1980s and developed ever since. Video 
cameras were used during several experiments, capturing wave 
run-ups or measuring sand bar movements [25,26]. Whereas 
traditional methods {in situ measurements, aerial and satellite 
imagery) are expensive and limited by weather conditions, video 
systems enable the researches to continuously monitor bea[:h 
processes in the long-term [27,28]. Especially since the technirjue 
has been automated and turned from analogue to digital, it has 
become a standard procedure when calibrating models for wave 
run-up and rip current occurrence [23,27,29]. 

Research described in this paper was conceived while 
examining possibilities to expand the beach inundation module 
oi SatBaltyk Operating System-Shores beyond the three test sites it 
currentiy includes. Preliminary analysis found that running 
XBeach model (ID variant) for an estimated 500 cross-sections 
of the coast (one for each kilometre of the coast) would be too 
time-consuming and therefore render the resulting forecasts 
useless. Moreover, it was designed specifically sandy coasts, which 
is not necessarily the case along the Baltic coast. Empirical 
equations of wave run-up on slopes were chosen as an alternative, 
for several reasons. Firstly, they are computationally undemand- 
ing, therefore creating a possibility to deliver forecasts without 
substantial delays. Secondly, they utilize the very same data as the 
system based on XBeach: wave parameters from WAM model, 
sea levels from M3D model and Maritime Offices' cross-sections 
of the shore. 

It excepted that the formulas, once calibrated, will be built into 
an operational, forecasting system during 2014—2015 in support of 
the SatBaltyk project [14]. Each 1 km cross-section wiU contain 
data extracted from the latest available high-resolution digital 
elevation models, created using lidar technology and provided by 
the responsible maritime offices. 

Additionally, Salhallyk project expects its system and models to 
be used, apart from delivering forecasts covering up to several 
days, also to archive the results and calculate hindcasts. They wiU 
allow us to conduct seasonal and long-term analysis. Therefore, 
alongside the calibration of the method, we include elements of 
such an investigation for two test sites, on which this experiment is 
based. 

Materials and Methods 

Wave run-up is defined as "the landward extent of wave uprush 
measured vertically from the stiU water level" [24]. The earhest 
formulation of this is by Hunt [30], who gives the following 
equation: 



where R is the wave run-up, Hg is the deep-water wave height and 
t,o is the deep-water surf similarity parameter. The latter is given 
as follows [31]: 



«0 = 



tan(a) 



(2) 



where a is the slope angle and Lq is the linear theory deep-water 
wave length: 



2n 



(3) 



where T is regular wave period. Hunt's formula can be used to 
calculate run-up from incidental regular waves. Later research 

reworked eq. (1) to include irregular waves and introduced 
additional empirical parameters based on field measurements and 
lab experiments [24]. Several very similar equations were 
developed, of which the most recent version from Mase [32] was 
chosen: 



J?2% = aH„ 




(4) 



R 

Ho 



(1) 



where R^x is wave run-up with a 2% probability of occurrence 

(i.e. 1 in 50 waves), H,^^„ is the significant deep-water wave height 
and L„p is the same as in eq. (3), except that T stands here for peak 
wave period (T^p). (The 2% run-up is commonly used in similar 
equations, though the reason for such a choice is not well 
explained. It seems to have a special statistical ridation to 
significant wave height and originate from Dutch experiments 
carried out in the 1930s, but the report where these results were 
published apparently has gone missing [17]). Additionally, two 
empirical parameters can be found in eq. (4): a and h. Mase 
suggested a = 1.86 and h = 0.7 1, basc-d on his laboratory tests [32]. 
Several field tests have been conducted by various authors, but 
they were mosdy located on high-wave-energy coasts of the United 
States, the Netherlands, Portugal, Spain or Australia 
[23,28,33,34]. One study for the Baltic sea could be identified, 
to our best knowledge, albeit it did not include any adjustmc-nts of 
parameters [35] . As a result, the formula had to be calibrated on 
the Polish coast to ensure- its accuracy. It should be also noted that 
other field-tested equations exist, developed by Hughes [22] and 
Stockdon [23]. Melby [24] tested those equations versus two 
incarnations of Hunt's formula using Stockdon's compilation of 
field experiments and found that Holman's [19] and Mase's 
equations had better statistical performance. Equation by Holman 
differs from Mase's only by the inclusion of a third parameter. 
Their performance, when calibrated, is virtually of the same 
quality, therefore Mase's variant was chosen as having less 
parameters to calibrate. This is even in spite of the fact that the 
main assumption behind the formula — a plain, impermeable 
slope — conforms more to flood protection structures than beaches, 
especially sandy ones. 

Pohsh Baltic Sea coast is 500 km long, including 72 km along 
the Hel Peninsula. Approximately 80% is a dune coast, whereas 
the remaining part forms a cliff coast [36]. Vast majority of the 
coastline has a beach, except for river mouths, harbours, military 
bases and the swampy shores of Puck Bay. Research was carried 
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Figure 1. Polish Baltic Sea coast. Location of the test sites in 
IVIi^dzyzdroje (1) and Dziwnowek (2) is marked. 
doi:10.1371/journal.pone.0105437.g001 

out in two locations in the western part of the coast, namely 
Dziwnowek and Migdzyzdroje (Fig. 1). 

Migdzyzdroje is a town on Wolin Island, while Dziwnowek is a 
village located 30 km east of the former. Both are very popular 
summer resorts along the banks of the Pomeranian Bight. The 
bight is a shallow basin, where depths are generally below 20 
meters. The coast in this area is oscillating, but in overall erosion is 
stronger than accumulation. It includes both cliffs and dunes, 
which are retreating, on average, by 1 0 cm per year [37]. 

The two selected cross-sections are located at about 414.2 and 
386.2 km of the Maritime Office's kUometrage. Both shores 
include moderately developed dunes that reach 5-7 m above sea 
level. There are no coastal protection structures in those particular 
locations apart from some timber groynes in Dziwnowek. They 
are also under continuous video surveillance. 

In Dziwnowek a video camera was mounted in 2009 in order to 
monitor a test site for the MICORE project [7]. Another camera 
was installed in Migdzyzdroje in 2013 to support model vahdation 
activities in the SatBaltyk project [14]. They are both placed at the 
tops of observation towers used formerly by the military and 
border guards. The view from the cameras is presented in Fig. 2 
(real-time images could also be found online [38]). The cameras 
were specifically adjusted for use as a close-range remote sensing 



device through a calibration processes described in literature 
[39,40]. 

The cross-sections under observation were surveyed with GPS- 
RTK devices while the cameras were operational, so that the exact 
location of each point along the profile could be unequivocally 
identified on all images from video monitoring. Measurements 
were conducted on July 11*, 2013 (Dziwnowek) and July 18*, 
2013 (Migdzyzdroje). In Dziwnowek the length of the profile, 
measured from the waterline at mean sea level to the base of the 
dune, is almost 33 m with an elevation difference of 2.32 m. For 
Mifdzyzdroje the values are 53 m and 1.88 m, respectively. 
Detailed data can be found in Table SIA in File SI. 

A major issue is calculating the slope angle properly. The 
method assumes a single value, whUe in reality the beaches are 
more complex. The profiles shown in Fig. 3 are relatively smooth 
and do not include a longshore bar that typically appears near the 
waterline after storms. Presence of a sand bar can skew the results, 
but it did not appear throughout the course of the experiment. 
Due to the limitations of the measurement method used, the cross- 
sections do not include most of the underwater part of the 
foreshore. Approaches to the slope angle problem differ largely. 
For example, Douglass [20] suggested skipping the slope 
parameter altogether, based on his field experiments on the US 
Adantic coast. On the opposite end, an 'effective slope' method 
was presented as the most accurate by Mayer and Kriebel [41], 
albeit it is most cumbersome to apply, as the slope angle is 
recalculated for every event between the wave breaking point and 
run-up on the slope (an iterative method). For the purpose of this 
study, the angle was calculated using linear regression based on all 
available measurements from the underwater part up to and 
including the dune base. Consequendy, the angle amounts to 3.9° 
in Dziwnowek and 2.0° in Mifdzyzdroje. Notwithstanding, large 
disproportion between beach height and width (1:12 and 1:23, 
respectively) means that the assumption of a plain slope is still 
fairly accurate, at least as long the beach is not reshaped by a 
storm. 

The final element of eq. (4) that had been acquired were data on 
waves, i.e. deep-water significant wave height and peak wave 
period. Since no actual measurements from buoys were available 
when the research described here was carried out, calculations 
from WAM wave model were used instead. This two-dimensional 
model of wave spectrum [42] was operationally implemented for 
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Figure 3. Cross-sections of two beaches in Mi^dzyzdroje and 
Dziwnowek. 

doi:1 0.1 371/joumal.pone.01 05437.g003 

the Baltic Sea by the Interdisciplinary Centre for Mathematical 
and Computational Modelling [43]. It provides a 84-hour 
forecasts every 12 hours [44] (ICM 2014). The latest forecast for 
each event included in the calibration process was used. ICM's 
model has a 8 km grid; central points of those grids utilized in this 
study are marked in Fig. 1. WAM point no. 12 supplied 
information for Migdzyzdroje, while point no. 19 for Dziwnowek. 

One additional factor had to be included, namely sea level, 
because eq. (4) calculates run-up from stUl water level. Measure- 
ments were obtained from the Institute of Meteorology and Water 
Management, which publishes real-time data from tide gauges on 
its website [45] . Levels from gauges in Dziwnow and Swinoujscie 
were used for the calibration process in Dziwnowek and 
Migdzyzdroje, respectively. The Baltic is a non-tidal basin, which 
also marks a difference from other research, which was done 
hitherto on tidal coasts. 

Wave run-up was extracted from two time periods during every 
day from June 20* up to December 8*, 2013. Such approach 
differs from other experiments [23], which analysed a few days or 
weeks in great detail. Here, almost half a year in less detail was 
chosen in order to ensure the method's validity over a longer 
period of time, given our ultimate goal of creating a forecasting 
tool. Only days with considerable wave heights or sea levels were 
analysed in one-hour periods, which was the temporal resolution 
of hydrologic data available. Such events occurred on 23*— 24* of 
September, 24" -26"' and 28" -29* of November, 2"" and 5*-7* 
of December. This last event was a fierce winter storm Xavier, 
which completely reshaped the beaches, creating a vast sandbar 
near the waterline and an indention towards the middle of the 
profde. In effect, calibration process was not continued after that 
date. During periods of calmer weather, wave run-ups from 
around 7 and 1 3 UTC were an analysed, with small alterations on 
some days due to availability of data. Thanks to this approach, 
small inundations of 0-20% would not completely dominate the 
calculations, even though they still represented the bulk of data 
acquired. 

For each event that was selected, up to 15 minutes of footage 
were observed in order to extract the highest wave run-up on all 
images. Since there would be about 100 individual wave run-ups 
during this period of time the value obtained from the images is 
likely to correspond to the 2% run-up modelled by eq. 4 [33,46]. 
Measurements of the highest run-up during in the 15-min set of 
images sampled at 1 Hz were made along the profde lines 
mentioned above. Distance between the point of the highest run- 
up and the 'zero' point could therefore be obtained from the 
image. Thanks to GPS RTK measurements at different points 



along the profile line, the distance in pixels could be easily 
transformed to distance in the terrain. As a result, a percentage 
value of beach inundation could be obtained. The 'zero' point (0% 
extreme wave run-up) was located at the point of the profile that 
equals the mean sea level, defined as —0.08 m in Kronstadt-86 
reference system or 500 cm at Polish gauge stations. On the other 
end of the scale, 100% run-up would mean that the extreme (2%) 
waves reached the dune base. The distance between those points 
constitutes the 'baseline' beach width. Therefore, modelled 
extreme beach inundation B could be written as: 



B=-^ — 100% 
h 



(5) 



where Rvv,, is the wave run-up from eq. (4), W is the sea level and h 
is the 'baseline' height of the beach. The inundation parameter 
could a have a negative value, since the sea level could be so low, 
that the waves reached below the mean sea level point. Therefore, 
a —10% value indicates that the beach is at the moment 10% 
longer than the 'baseline' width. On the other hand, maximum 
values of 5 were capped at 100%, since any higher value means 
that the waves reach outside the beach. 

A total of 667 events were analysed: 367 in Mi(;dzyzdroje and 
300 in Dziwnowek (the latter number is lower due to camera 
malfunctions). During those events, water levels varied from 
426 cm to 596 cm; significant wave heights reached from 0.01 to 
3.99 m and peak wave periods from 2.0 to 10.2 seconds. Large 
variety of values indicate that both storm surges and falls were 
included, as well as period of calm weather. The lowest measured 
beach inundation amounted to —22.5% in Migdzyzdroje and — 
7.3% in Dziwnowek, while during the Xavier storm waves 
swashed the dunes in both locations (100% run-up). Histograms 
are presented in Fig. 4. Underlying data (except for sea levels) can 
be found in Table SIB in File SI. 

Having obtained actual wave run-ups, empirical parameters a 
and b where adjusted, so that the linear regression line plotted 
using the data for each locality would pass through points of equal 
value of both actual and modelled run-up, i.e. constitute a 1: 1 hue. 
The correlation between the modelled and actual data was not 
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Figure 4. Histograms of parameters included in the study. 
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maximized liere, but it decreased the formula's bias in estimating 
high run-ups. 

A permit was required in order to conduct GPS measurements 
on the dunes. It was obtained from the Maritime Office in 
Szczecin. No other permissions were necessary. Video cameras 
used in this study were placed on the university's (Mig dzyzdroje) or 
private (Dziwnowek) parcels by kind consent of their respective 
owners. 

Results 

Results of the calibration process are presented in Fig. 5. 
Empirical parameters a = 1 .29 and 6 = 0.72 where used, so that eq. 
(4) becomes: 



i«2% = 1.297/,,, 




0.72 



(6) 



In general, the data concentrate along the regression line, with 
very few oudying elements. Root-mean-square error for data from 
Migdzyzdroje is 6. 1 percentage points, or 3.2 m of beach width. In 
Dziwnowek, the error is slightly larger, 6.5 percentage points, but 
because the beach is shorter, this value corresponds to 2.2 m of its 
width. The coefficient of determination is high and amounts to 
0.94 in Migdzyzdroje and 0.92 in Dziwnowek. 
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Figure 5. Measured and modelled beach Inundation for 
Mi^dzyzdroje and Dziwnowek. 
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In Migdzyzdroje, the formula proved most accurate for small 
run-ups, up to 30%, while more energetic events where slighdy 
underestimated. In Dziwnowek, it is the opposite. However, these 
are only slight deviations from the line of equality. On most 
occasions. Unfortunately, only a few cases of 60-80% run-ups 
where recorded, causing some uncertainty about the method's 
reliability for storm events, even though ca. 90% beach inundation 
was estimated correctly. It generally properly predicted waves 
reaching the dune base. 

However, the impression of very good quality of the model is 
slightly misleading. An important factor regulating how far will the 
wave swash into the beach is the sea level. It contributed, on 
average, to about 22% of the actual beach inundation in 
Mifdzyzdroje and 32% in Dziwnowek. Since the data from tide 
gauges are fairly accurate (even though there is some considerable 
distance from the gauge station to Mijsdzyzdroje), performance of 
'bare' wave run-ups was calculated. Normalized RMS error 
amounted to 0.27 in Mifdzyzdroje (i. e. average error in 
estimating wave run-up was 27%) and 0.29 in Dziwnowek. 
Coefficient of determination was 0.85 and 0.81, respectively. This 
is a lower value than that recorded by the beach inundation 
indicator, but stiU a decent result. 

Having obtained optimal parameters for the formula, wave run- 
up was recalculated using hourly data from WAM and M3D 
models. The latter is a hydrodynamic model providing, among 
other parameters, water levels for the Polish coast [47,48]. Fig. 6 
presents a histogram of inundation from hourly data and its 
change over the course of 2013 (daily averages are presented for 
clarity of the graph). Run-ups below 30% constituted four-fifths of 
all events. The values for the beach in Dziwnowek were generally 
higher than in Migdzyzdroje (20% against 12%). The reasons 
could be easily identified: not only Dziwnowek has a shorter and 
steeper beach, but also had higher waves {H,„„ = 0.53 m against 
0.47 m) and water levels (504 cm against 500 cm) during that 
period. The beach in Migdzyzdroje was more inundated than 
Dziwnowek' s only 5% of the time. 

In Fig. 6 both calmer and more dynamic periods are visible. 
Especially standing out is the Xavier storm surge in early 
December. Only one event in February comes close. November 
and December obseived the highest monthly means, with April 
and October on the opposite. In Dziwnowek, 18 days had average 
with inundation above 50%, whereas in Migdzyzdroje the 
corresponding figure is 9 days. 2013 was not a particularly 
eventful year in terms of storms in general. Even Xavier, though 
generated strong winds and caused damages in some Western 
countries was at the end rather unimpressive in terms of coastal 
ffooding. It could be estimated that it impacted the dunes in 
Dziwnowek for around one and a half day, and even less than that 
in Migdzyzdroje. Naturally, the results for parts of the year could 
be not exactly accurate due to changes in morphology of the 
beaches, particularly after the December surges, where sandbars 
near the waterline appeared. 

Discussion and Conclusions 

The empirical parameters' values from eq. 6 do not differ much 
from other authors' findings. Parameter b, estimated at 0.72, is 
almost the same as Mase's lab-tests result of 0.7 1 [32] and Melby's 
approximation of 0.70 based on Stockdon's dataset of nine 
different field experiments [23,24]. Such a small difference 
between different studies could indicate that the value of around 
0.7 has a universal application. On the other hand, a parameter's 
value (1.29) is significantly lower than Mase's 1.86 [32]. Other 
authors' estimates of this parameter ranged from 1.49 to 1.87 [24]. 
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Figure 6. Extreme 2% beach inundation estimates for iVIif dzyzdroje and Dziwnowek, hourly data from January 1 to December 

31 2013. Daily averages (A) and histogram (B). 
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It is likely caused by a different research scope: Mase and others 
tested their methods using very large slope angles that are typical 
for hydrotechnic structures such as dykes. In contrast, Melby's 
estimate fi'om field data is 1 . 1 0 [24] . It is likely caused by different 
environmental conditions, since the beaches used in that study 
were located along tidal coasts with high wave energies. Another 
source of the difference may be the permeability of the beaches, 
though authors working on the empirical equations generally did 
not give much consideration to this aspect apart from assuming 
that the slopes they tested were impermeable or permeable, 
therefore is hard to quantify the impact of this factor. Overall, 
bigger permeability lowers run-up [49] and the Polish beaches are 
definitely more permeable than those used in other studies. 

Our results could be compared with other findings using 
Melby's Summary Performance Score. This statistic is an average 
of normalized RMS error, normalized bias and scatter index with 
an 0 to 1 scale, where 1 means no error [24] . Only wave run-up 
from stiU water level is taken into account here. In our case, the 
Performance score amounted to 0.78 in both locations. Applica- 
tion of Mase's equation adjusted by Melby to Stockdon's field data 
gave a result from 0.62 to 0.90 depending on the test site. Four out 
of nine beach had better score than 0.8, with an average of 0.84 
[24]. Melby's findings were more accurate than our results, though 
it should be remembered that they were based on a much smaller 
number of run-up events (from 14 to 138, depending on the test 
site). Additional error could have come from the limitations of our 
dataset: the use of modelled wave data and inaccuracies in 
registering wave run-ups from images. Nevertheless, it can be 
concluded that the results of the calibration process were in line 
with other authors' findings. 

Notwithstanding, some shortcomings of the method could be 
identified in the context of application for the whole Polish coast in 
SatBaUyk Operational System-Shores. First of aU, the two test sites 



are both located on the western part of the coast, therefore not 
necessarily accurate enough for the middle and eastern section. 
Currently there is no possibility of validating the model there, since 
no suitable, scientific cameras are installed along the coast. What is 
more, the method has its potential inaccuracies originating in the 
heavy reliance on other models' output: modelled wave data are 
used for both calibration and operational forecasting. AU this 
means that the method described here, though accurate for the test 
sites, should be subjected to further validation in other areas, 
before implemented to the entire Polish coast in an operational 
system. 

When completed, the forecasting system could have several 
applications. For example, temporary constructions that are 
placed on the beach in the summertime could be removed in 
time to protect them from the effects of high wave run-up. 
Moreover, run-up predictions close to 100% would indicate 
potential dune/cliff erosion, an information useful for scientists 
interested in the subject of coasdine change. As could be seen in 
Fig. 6 the method is also useful for analysing storm surges' 
development over time, as well as seasonal analysis. The topic wiU 
be investigated further in the future. More work is also expected 
on the calibration, too, pending availability of new data. 

To sum up, the study proved that empirical wave run-up 
formulas could be applied to the Baltic Sea after calibration to 
local conditions. The method is accurate enough to be used in 
other work, especially lays foundations for an operational, 
forecasting system (which wiU naturally depend also on the quality 
of operational hydrologic models delivering data for calculating 
beach inundation). The errors recorded were not substantial 
despite applying a formula designed for very different conditions. 
This proves that this simple method, if properly used, can also be a 
valuable addition to more complex research. 
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